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Abstract. 

We present electronic transport measurements through thiolated Cgo molecules 
in liquid environment. The molecules were placed within a mechanically controllable 
break junction using a single anchoring group per molecule. When varying the electrode 
separation of the Cgo-modified junctions, we observed a peak in the conductance traces. 
The shape of the curves is strongly influenced by the environment of the junction as 
shown by measurements in two distinct solvents. In the framework of a simple resonant 
tunneling model, we can extract the electronic tunneling rates governing the transport 
properties of the junctions. 
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1. Introduction 

How can we experimentally determine the relevant energy levels and the tunneling rates 
involved in electronic transport through single molecule devices? Such questions are of 
fundamental interest for the development of molecular electronics. The conductance of 
single molecules has been investigated in several experiments where phenomena such 
as Coulomb blockade and Kondo effect [1, 2], negative differential resistance [3] and 
logic gates [4] were pointed out. Different techniques were used to address single 
or a few molecules. Fixed contact arrangements such as nanopore systems [5] and 
nanogaps fabricated by shadow evaporation or electromigration [6] have been realized. 
Alternatively, scanning tunneling microscope (STM) [7, 8] and break junction techniques 
[5, 9, 10] offer tunable contacts to characterize molecules. Break junctions offer an 
additional stability when compared to the STM, but the implementation of a third 
electrode remains delicate although recent work shows it can be done [11]. Adding a 
chemically controllable environment in break junctions experiments would provide an 
additional control of the anchoring of the molecule to the metallic constriction and allow 
in situ electrochemistry experiments. The excellent efficiency of a liquid gate bas been 
demonstrated previously on carbon nanotube field-effect transistors [12], and recently 
this effect has been studied for organic molecules using a STM [13, 14]. 

Buckminsterfullerenes and in particular Ceo an d its derivatives have attracted much 
attention [15] and, due to their particular electronic properties, represent ideal model 
systems for molecular electronics [16, 17, 18]. In this Letter, we present electronic 
transport measurements through thiolated C 60 derivatives in a two probe configuration 
provided by the tips of a break junction (Fig. 1). The single functional group serves 
as an anchor and, in the case where a single molecule junction is realized, permits to 
control the coupling of the Cgo to a Au electrode by mechanically adjusting the inter- 
electrode spacing d. By performing the transport measurements in two distinct solvents, 
dimethyl sulfoxide (DMSO) and toluene, we could observe the strong influence of the 
environment on the molecular junctions. 

2. Experimental section 

The synthesis of the fullerene derivative was carried out following well-established 
protocols [19]. A thiolated C4 linker ensured that the C@o molecule would bind to 
the gold electrodes. To contact the molecules and perform transport measurements, 
we used a mechanically controllable break junction setup [20]. A gold nanobridge was 
lithographically patterned on a phosphor bronze substrate. By bending the substrate 
via a vertically moving rod, the gold bridge is elongated and finally breaks. The two 
resulting ends serve as electrodes for contacting the molecules. The setup includes a cell 
that allows working in different liquid environments and has been described previously 
[21]. 

We first characterized the gold junctions via measurements of the conductance 
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G(d) versus distance d in air and in the solvents used subsequently for the C 60 
molecules, toluene and dimethyl sulfoxide (DMSO). Several sets of consecutive open- 
close cycles were recorded, with G ranging typically between 0.001 Go and a few Go, 
where Go = 2e 2 /h is the conductance quantum. To investigate the transport through Cqo 
molecules, a 0.1 mM Ceo solution was added to the liquid cell while the junction was kept 
close. The junction was then opened widely 5 nm) to favor the self-assembly of the 
molecules. We also added molecules while keeping the junction open but no significant 
differences in the transport measurements were observed. Sets of 10 — 30 G(d) curves 
were recorded at time intervals of 30 — 45min. During these cycles, the junction was 
never fully closed, keeping its conductance below 0.1 Go in order to limit mechanical 
rearrangements of the Au tips. G(d)-values obtained while closing the junction are 
shown in Fig. 2 for DMSO (a) and toluene (b) at a bias voltage of V& = 0.2 V. 

The open symbol curves show measurements in the pure solvent. The solid 
curves are five successive raw measurements performed after adding the C 60 solution 
to the junction. We also attempted to measure current-voltage (I — V) characteristics. 
However, these curves were not reproducible, displaying rather large uncontrolled 
hysteresis, which we assign to the dynamics of the liquid environment and the relatively 
weak tethering of the C 60 with one single linker group only. 

While the measurements of G(d) in the pure solvents show a clean exponential 
behavior [21], the curves for the Ceo modified junctions exhibit a more complex structure. 
In previous break-j unction work with molecules, the occurrence of a plateau in G(d) has 
been taken as the signature for the 'locking-in' of a molecule [10] . We observe a plateau in 
about 50 % of the measured closing cycles. In addition to this observation, a pronounced 
peak can appear (arrows in Fig. 2). Though, the peak appears less frequently in different 
samples, reproducible measurements (visible in each closing curve within a cycle of 10 
curves at least) could be acquired for both C 60 in DMSO and toluene. In DMSO, the 
peak corresponds to a conductance maximum of ~ 0.012 Go and is followed by a plateau- 
like feature. When further closing the junction the conductance rises more rapidly, 
although less sharply than in the pure solvent. The shape of the conductance curves in 
toluene (Fig. 2b) is very different, showing more noise and a fully developed peak. After 
a small pre-peak (not present in all curves), the conductance rises to a maximum value 
of ~ 0.025Go at the peak and decays again to a very low conductance. While further 
closing the junction the conductance increases sharply, with a slope similar to that of 
the conductance in the pure solvent. The environment in which the junction is studied 
plays clearly an important role here. We note in particular that DMSO is a rather poor 
solvent for neutral fullerenes whereas toluene is a good one (see e.g. [22]). 

To calibrate each junction, we measured the tunneling conductance between the 
two Au electrodes in the pure solvent. Let us denote the gap distance between the 
two Au electrodes by d gap . It can be expressed as d gap = rs, where s denotes the 
vertical displacement of the pushing rod and r the reduction factor which depends on 
the geometry of the junction [20] . It is the reduction factor r which needs to be calibrated 
for each junction. The tunneling conductance G(s) measured in a pure solvent was fitted 
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according to the expression G(s) oc exp[— nrs] with k = 2y / 2m0//i, where <f> (and k) 
represents the height of a square tunneling barrier of width d gap and m is the electron 
mass. The fits were adjusted by tuning r in order to obtain values for <fi in agreement 
with our previous work [21] in which the apparent workfunctions in different solvents 
were calibrated with respect to measurements in vacuum. The reference inverse decay 
lengths k, valid for the pure solvents, are K to iuene = 0.85 A -1 and k,dmso — 0.63 A -1 . 
This procedure provides the reduction factor r for each junction measured. [23]. 



3. Discussion 



The particularity of our system lies in the single anchor group used to tether the molecule 
within the break junction. This configuration makes it possible to tune the coupling 
of the molecule to the electrode: while one tunneling rate (Yi/h) is fixed, the second 
(Y 2 /K) varies exponentially with the gap between the electrode (electrode to the right 
in Fig. 1) and the C 60 molecule. The metal-molecule-metal junction forms a double- 
barrier junction. In this picture, the observed conductance peaks can be explained 
within a coherent resonant tunneling model, provided that one of the energy levels of 
the molecule lies not too far away from the Fermi level of the electrodes (Fig. 3a). 

Resonant tunneling through a single level e results in a peak in the linear 
conductance G when the level aligns with the Fermi level Ep of the leads, i.e. at 
e = 0. The occurrence of a peak in G(e) is a well known fact [24]. However, a peak can 
also arise in G(Y 2 ), even if e ^ 0. Resonant tunneling is described by the Breit-Wigner 
equation: 

i*e 2 + (r 1 + r 2 ) 2 /4 1 1 

In our geometry, Ti is fixed and defined by the molecular tether holding the 
molecule to the left electrode (Fig. 1). T 2 depends exponentially on the distance d 
between the C 60 molecule and the right electrode with k being the inverse decay length, 
i.e. Y 2 (d) = Y^ ■ e~ Kd . G(Y 2 ) in Eq. (1) as a function of T 2 has a maximum G max 



at = A/4e 2 + rf. For e « Ti we would have = I\ and consequently a peak 
conductance of G max — Gq. In the experiment, the peak conductance is however much 
smaller than Go, so that Ti <C e. It follows then that — 2e and G max — Gq-Y\/e. A 
respective plot in this limit (ri <C e) is shown in Fig. 3b. Note, that we define e > 0, 
but the sign of e is not determined, i.e. a sign change in e does not affect the discussion. 

Using the Breit-Wigner equation we can obtain the ratio Y\/e from the measured 
peak height without fitting. The measurements yield Yi/e ~ 1.1 • 10~ 2 in DMSO and 
Yi/e ~ 2.7- 10~ 2 in toluene, more than twice the DMSO value. Because we expect Ti to 
be the same for both solvents (after all it is the same linker), this suggests that the level 
is lying closer to the Fermi level of the electrodes in toluene as compared to DMSO. 

That we can measure directly the ratio Y\/e is a very nice fact. But we would like 
to obtain both e and Yi. This is at first sight not possible, because Eq. (1) is scale- 
invariant. Changing Ti, Y 2 , and e by the same factor leaves the curve invariant. To our 
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rescue, we emphasize that we did not measure the linear-response conductance, but G 
at a bias of V& = 0.2 eV, much larger than ksT at room temperature. The finite bias 
introduces an energy scale which allows a considerable narrowing down of the range of 
possible values. 

Integrating Eq. (1) over the bias window, given by the applied voltage Vj,, yields: 



= 4e 2 7i 72 (d) 
h 71 + 72(d) 



2e \ f l + 2e 

arctan — + arctan 



7i + 72(d)/ \7i + 12(d) 



(2) 



All parameters are now scaled to i.e. 7$ = Ti/eVt, 72(d) = 72 • e~ , and e = e/eVj,- 
There are two regimes: a) e < Vb/2 and b) e ^> V&/2, the latter corresponds to the Breit- 
Wigner equation. In both cases a peak develops in £(72), but there is a pronounced 
crossover, which sets in sharply at e := e/eV b = 0.5. This crossover is shown in Fig. 3c 
for a selected value of 71 = 0.02. In the following, we will try to fit our measurements 
to Eq. (2) in both the low and large e limit. There are four fitting parameters: e, T 1: 
r^, and k. Note, that we always use the ansatz T 2 = • e~ Kd and that we define the 
zero of the distance axis d to match with the peak position in the measurement. This 
is possible, as we have no means to determine the 'true' zero, i.e. the point of contact 
of the right electrode with the Ceo molecule. 

In the first step, we obtain k by focussing on large d where the junction is widely 
open and the conductance well below the conductance at the peak. In this simple 
tunneling regime, G oc Y 2 oc e' Kd . From fits to the averaged conductance curves 
of the C 60 -modified junctions (Fig. 2 insets), we extract the inverse decay lengths 
k-dmso = 0.43 ± 0.01 A -1 and K to i uene = 0.37 ± 0.05 A -1 . We emphasize that the 
value obtained for DMSO is more reliable than the one for toluene, since the fit could 
be performed over a more extended range of conductance. This is due to the presence 
of a lower conductance peak in toluene (see Fig. 2b). The inverse decay lengths k are 
slightly suppressed as compared to the values measured in the solvent alone [21]. 

We next compare our data to conductance curves calculated in the limit e — 0. We 
fix the k values to the ones determined for large d and try to obtain a good match 
with the experiment by tuning both Ti and I^. The dashed lines in Fig. 4 show 
curves for ri )DM5C , = 0.92 meV and r* iDMSO = 2.7 meV, and T 1Muene = 2.2 meV, 
^2 toluene = 3.9 me V. In both cases, the adjusted curves present a substantially broader 
peak than the data, showing that the limit e = is inappropriate. 

We now relax the condition e = and try to find the smallest possible e for which 
a good match between the calculated and measured curve is obtained, again using the 
k values determined for a large d. We obtain Sdmso — 0.20 eV, T 1DM so — 4.3 meV, 
t *2,dmso — 0.35 eV, and e to i uene ~ 0.12eV, Ti Muene ~ 5.5 meV, T* 2 toluene ~ 74meV. 
Care has to be taken at this point, because Etoiuene — 0.12 eV is very close to the 
transition region between the Breit-Wigner and the e = case where the theoretical 
G(d) curves are quite sensitive to the temperature. When taking the finite temperature 
into account, e to iuene slightly increases, so that e > 0.15 — 0.2 eV is a good approximation 
for both DMSO and toluene. 
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If we allow for larger values of e we cross over to the Breit-Wigner equation. In 
this scale-invariant limit, we obtain convincing fits for both DMSO and toluene with the 
parameters Y\/e = 0.011 (DMSO) and 0.027 (toluene), and T 2 /e = 2 as required (solid 
lines in Fig. 4). Hence, at this stage we can say with confidence that e > 0.15 — 0.2 eV 
for both DMSO and toluene. 

As emphasized before, the Breit-Wigner formula cannot provide an upper bound 
for e. However, we can - on a physical ground - provide an upper bound on T 2 and 
therefore also on e. The latter is obtained through the relation = 2e. The peak 
in G(r 2 ) appears at relatively large T 2 values, but T 2 cannot grow to too high values, 
because otherwise higher molecular states of the Ceo molecule will strongly hybridize 
with the states in the Au electrode and the model will break down. Moreover, a large Y 2 
at short distances corresponds to a large binding force of order kT 2 between the molecule 
and the Au electrode. In fact, we think that the crossover from the peak to the plateau 
in DMSO is determined by this consideration. Furthermore, the fully developed peak 
in toluene suggests that T 2 at the position of the peak (T 2 ) is smaller in toluene than 
in DMSO. 

In order to fix the upper bound of T 2 we need to discuss the typical energy scales 
of the orbitals in C 60 [25]. In the inset of Fig. 4 we show the energy levels of C 60 and 
C 60 (adapted from Green et al. [26]). The HOMO-LUMO gap of C 60 determined in 
solution is around 2.3 eV [27] whereas local density approximation calculations provide 
a value of 1.9 eV for the free Ceo molecule (see e.g. [15] and refs. therein). Experiments 
on Ceo monolayers assembled on metallic surfaces showed values ranging from 1.5 eV to 
3.0eV [28, 29, 30, 31]. 

Due to its high electronegativity, a C 60 molecule in contact with a metal tends to 
gain charge from the metal. The ionization potential of Au is 5.3 eV and the electron 
affinity of C 60 equals 2.8 eV. The latter includes the Coulomb energy of the singly 
charged Ceo amounting to ~ 3 eV. If we assume that the Coulomb energy is strongly 
screened in the gap between the Au electrodes, C 60 ^Cg is favorable. Photoemission 
studies of Ceo on a Au surface yield values for the electron charge on Ceo of 0.8 [32] and 
1.0 [33]. Hence, it is very likely that Ceo is singly charged also in our work. On acquiring 
an electron the threefold degenerate LUMO state of Ceo splits by a very small amount of 
8 meV (not shown and also not resolvable in our experiment) and the chemical potential 
E F is expected to align closely with this orbital (dashed line). In fact, E F may lie above 
(e > 0) or below (e < 0), depending on the exact charge state. 

Taking the HOMO-LUMO gap of the C 6 o and the C 6 o anion as the proper 
energy scales, T 2 < 1...2.3eV to avoid hybridization with higher lying orbitals. If 
we impose that this condition occurs in DMSO at the transition from the peak to the 
observed plateau at d — —2 A (see arrow in Fig. 4), we obtain £dmso < 0.2 . . . 0.5 eV. 
Together with the previous consideration, e is determined quite accurately to 
0.2 < Edmso < 0.5 eV. Similarly, if we pretend to describe the full peak in toluene 
up to d = —7 A (arrow), we obtain e to iuene < 0.04 . . . 0.09 eV. These rather small values 
for e are inconsistent with the previous consideration. In order to obtain an upper 
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bound for e to i uene of at least 0.2 eV, one would have to set the cut-off to as large a 
value as T 2 < 5 eV. This disagreement shows that the full decay of G(d) in toluene to 
approximately zero conductance beyond the peak must have another origin. This will 
be addressed in the following. 

We have assumed here that we are allowed the use the resonant tunneling model 
with an unperturbed level beyond the peak. However, a full understanding of the small 
d regime would require the understanding of local rearrangement effects taking place 
when the electrodes and the molecules are reaching mechanical contact. In the case 
of toluene we observe that the conductance drops to a negligibly small value after the 
peak and then shows a pronounced increase, similarly to the behavior of the junction in 
the pure solvent (Fig. 2b). We can presume that, upon closing of the junction to small 
values of d, the molecules are pushed away from the inter-electrode space, resulting in 
conductance properties dominated by Au-Au tunneling. The shoulder observed in the 
case of DMSO (Fig. 2a) lets envision a different scenario in which the molecules remain 
within the junction and therefore continue to substantially affect the conductance of 
the junction at very short inter-electrode separations. Additionally, large conductance 
fluctuations can be observed in toluene (see Fig. 2b) while the curves are much smoother 
in DMSO. We attribute these differences to the solubility of the fullerene derivative in the 
solvent considered. Whereas toluene is a good solvent, the fullerenes tend to aggregate 
with time when suspended in DMSO. In our view, the large conductance fluctuations 
observed in toluene reflect the tendency of the fullerene molecules to be re-solved if the 
cavity is getting to narrow. In contrast, the molecules rather stay within the junction 
if DMSO is used due to the lower solvation energy. 

Whether a single molecule picture really holds in such a system is an important 
question. The very good match between theory and experiment supports the single 
molecule approach. The analysis provides more than a qualitative agreement with the 
data and allows the extraction of narrowed range estimates for the LUMO position and 
the tunneling rate Ti of the molecular bridge between the Au electrode and the Cgo 
molecule. The results are: e = 0.2 . . . 0.5 eV, and Fje =11... 2.7 • 10" 2 . 

The above discussion is based on a non-interacting model, where the interaction 
energy has been absorbed in the level energy e. When the coupling r 2 to the molecule 
changes, the average charge on the molecule is modified, thereby changing the interaction 
contribution to e. However, we have numerically checked (within a self-consistent 
mean field theory [34]) that this effect does not appreciably change the shape of the 
conductance versus distance curves, but only slightly scales the parameters e, F 1 and 

1 2- 

The presence of a much lower peak in some raw conductance curves in toluene 
(about 12 A prior to the main peak, Fig. 2b) might be a signature for a weaker tunneling 
process involving Ceo molecules on both electrodes. If running experiments confirm this 
feature, an interpretation in terms of a triple-barrier system could further support this 
approach. 



Resonant tunneling through a Cqq molecular junction in liquid environment 



8 



4. Conclusions 

We have performed transport measurements through Ceo molecules tethered by a single 
anchor group in Au break junctions operated in a liquid environment. The signature of 
the presence of the molecules within the junction appeared as a peak in conductance 
versus distance curves. The shape of this peak was strongly influenced by the solvent 
in which the junction was operated, thereby showing the importance of a proper 
environmental control in nanoscale junctions. 

The data can be understood in a resonant tunneling picture through a single level 
located close to the Fermi energy Ep- The particular geometry of the experiment 
combined with the environmental control allowed us to provide numbers for the 
electronic tunneling rates between a gold electrode and a fullerene derivative. By 
systematically varying the linker group, our approach will allow to better understand the 
electronic coupling between molecules and electrodes in molecular devices. Moreover, 
gating will allow the resonant level to be moved even closer to the Fermi level. An 
increase of the conductance up to the quantum conductance value can then be envisaged. 
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Figure 1. Schematic representation of a break junction with a thiolated Ceo molecule 
anchored to the left electrode. The distance d between the molecule and the right 
electrode can be adjusted by opening and closing the junction. 
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Figure 2. Conductance G versus distance d curves in DMSO (a) and toluene (b) 
measured at a bias voltage of 0.2 V while closing the junction. The open symbol 
curves were measured in the pure solvent. The solid lines are raw data for five successive 
measurements after adding the Cgo solution. The insets show averaged (over 10 curves) 
conductance data in a log-lin representation. The value for the inverse decay rate n 
was obtained from a linear fit (solid lines) to G/Gq in the low conductance regime: 
k-dmso = 0.43A -1 and Ktoiuene = 0.37A" 1 . Note, that the zero in the distance scale 
d was chosen to match with the position of the peak. 
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Figure 3. (a) Schematics of the energetics of a double-barrier junction with a 
single molecular level positioned at e relative to the Fermi energy Ep of the metallic 
reservoirs, (b) Linear-response conductance G as a function of T2 for the double-barrier 
shown in (a), assuming coherent tunneling and Ti <C s. (c) Same as (b), but taking 
into account the finite bias voltage V& 3> fc&T. This illustration has been obtained for 
a specific value of T% = 0.02 Vb- A pronounced crossover is observed at e = eVb/2. 
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Figure 4. Averaged conductance versus distance curves for DMSO (□) and toluene 
(o). The solid lines are fits to the data using Eq. (1) with Ti/s =1.1- iCT 2 (2.7 • 1(T 2 ) 
and k = 0.43A- 1 (0.37A" 1 ) for DMSO and (toluene). The dashed lines were calculated 
from Eq. (2) in the limit e = 0: T ± = 0.92 meV (2.2 meV), T* 2 = 2.7 meV (3.9 meV), 
and K = 0.43 A" 1 (0.37 A^ 1 ) for DMSO and (toluene). 

Inset: Schematic energy diagram of the metal-molecule-metal junction showing the 
energy level of Cgo and Cg . Due to its large electronegativity, Ceo tends to gain 
charge from a nearby metallic surface as illustrated by the position of its LUMO level, 
lying below the Fermi energy of the Au electrodes. 



